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The Standard Model augmented by the presence of gauge-singlet right-handed neutrinos proves 
to be an ideal scenario for accommodating nonzero neutrino masses. Among the new parameters of 
this "New Standard Model" are right-handed neutrino Majorana masses M. Theoretical prejudice 
points to M much larger than the electroweak symmetry breaking scale, but it has recently been 
emphasized that all M values are technically natural and should be explored. Indeed, M around 
1 — 10 eV can accommodate an elegant oscillation solution to the LSND anomaly, while other 
M values lead to several observable consequences. We consider the phenomenology of low energy 
(M < 1 keV) seesaw scenarios. By exploring such a framework with three right-handed neutrinos, 
we can consistently fit all oscillation data - including those from LSND - while partially addressing 
several astrophysical puzzles, including anomalous pulsar kicks, heavy element nucleosynthesis in 
supernovae, and the existence of warm dark matter. Furthermore, low-energy seesaws - regardless of 
their relation to the LSND anomaly - can also be tested by future tritium beta-decay experiments, 
neutrinoless double-beta decay searches, and other observables. We estimate the sensitivity of such 
j^Tj' probes to M. 
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I. INTRODUCTION 



Neutrinos have provided our first glimpse into physics beyond the Standard Model of elementary particles (SM). 
Neutrino oscillation experiments have given us unambiguous evidence that the three active neutrinos have mass and 
J> mix (see [l| for recent reviews of neutrino theory and phenomenology). As most important discoveries, the results of 
these experiments have raised more questions than they answered. Even from our limited knowledge of the neutrinos, 
it is clear that their properties, including sub-eV masses and large mixing, are quite different from their charged 
fermion counterparts. The true explanation of this puzzling behavior likely rests on the fact that neutrinos are the 
only known fundamental neutral fermions, but the exact reason behind this is still open to speculation. Neutrino 
masses, as deduced from oscillation experiments, hint at the existence of right-handed neutrino states, lepton number 
■ violation, new sources of CP violation, as well as a new energy scale. 

The seesaw mechanism Q is an appealing way to generate the observed neutrino masses and lepton mixing matrix. 
The idea is simple. Add an arbitrary number of singlet fermion states to the SM matter content. The triviality of 
their quantum numbers allows them to have Majorana masses of magnitude M, as well as couple to the SU(2)l 
Oh! lepton and Higgs doublets. The latter vertices become Dirac mass terms of magnitude fi after electroweak symmetry 
breaking. The standard theoretical prejudice is that the Dirac masses are of order the charged fermion masses, while 
the Majorana masses are at some very high energy scale, M 3> 100 GeV. If this is indeed the case, the resulting 
propagating neutrino degrees of freedom separate into two quasi-decoupled groups: mostly active states with very 
small masses m ~ ^ 2 /M suppressed by the new physics scale, and mostly sterile states with very large masses M. 
In this scenario, the mostly right-handed states are not directly observable. Indeed, it is possible that if such a 
high-energy seesaw is realized in Nature, its only observable consequence is that the mostly active neutrinos have 
mass and mix. 

Of course, there is no direct evidence that M - which we refer to as the seesaw scale - is large. Large M values are 
attractive for several reasons, including the fact that one may relate M to the grand unified scale. On the other hand, 
all M values are technically natural, given that when M vanishes the global symmetry structure of the Lagrangian is 
enhanced: U(1)b-l is a symmetry of the Lagrangian if M = 0, so that M is often referred to as the lepton number 
breaking scale. This point was recently emphasized in Recent analyses have also revealed that there are several 
low-energy choices for the seesaw energy scale that allow one to address outstanding problems in particle physics and 
astrophysics. The main reason for this is that, unlike in the high-energy seesaw, in a low-energy seesaw the mostly 
right-handed states do not decouple but remain as kinematically accessible sterile neutrinos. 

The data reported by the LSND short baseline neutrino oscillation experiment |4j can be explained by postulating 
the existence of light (m ~ 1 — 10 eV) sterile neutrino states. This result is currently being tested by the Fermilab 
MiniBooNE experiment [BJ and, if confirmed, will lead the community to seriously contemplate the existence of light, 
SM singlet fermions. It was pointed out in (see also |||) that if M ~ 1 — 10 eV, the right-handed seesaw neutrinos 
can easily play the role of the LSND sterile neutrinos. There is also evidence for mixed sterile neutrinos at other 
energy scales: eV sterile neutrinos aid in heavy element nucleosynthesis in supernovae, keV sterile neutrinos can help 
explain the peculiar velocity of pulsars, and remain viable warm dark matter candidates. In the past several months, 
it has been shown that the seesaw right-handed neutrinos may play the role of all these astrophysically/cosmologically 
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inspired sterile neutrinos 0, II] • 

In this paper, we consider the phenomenology of low-energy (M < 1 keV) seesaw scenarios, extending the analysis 
performed in : 3j in several ways. In Sec. |n] we review the generation of neutrino mass via the seesaw mechanism 
and apply it to relatively light right-handed neutrino states. We pay special attention to the most general active- 
sterile seesaw mixing matrix, whose parameters are the main object of our study. In Sec. IIIII we review the several 
different "evidences" for sterile neutrinos, and discuss whether these can be "fit" by the low-energy seesaw. We 
concentrate on exploring solutions that can accommodate at the same time the LSND anomaly and the astrophysical 
processes outlined above, but also discuss different combinations of the seesaw parameters capable of explaining only 
the astrophysics-related observables. In Sec. IIVI we examine other experimental probes that can be used to explore 
low-energy seesaws - regardless of their relationship to the LSND anomaly, pulsar kicks, and warm dark matter. We 
concentrate on the prospects of current/future tritium beta-decay and neutrinoless double-beta decay experiments. 
We conclude in Sec. E] by summarizing our results, commenting on the plausibility of this scenario, and offering a 
general outlook for the future. 



II. THE SEESAW MECHANISM AND ELECTRON VOLT NEUTRINO MASSES: PRELIMINARIES 

In order to account for nonzero neutrino masses, we add to the SM particle content three SU (3) c x SU (2)^ x U (l)y 
gauge singlet fermion states Ni, conventionally referred to as right-handed neutrinos. While sterile under SM gauge 
interactions, right-handed neutrinos may still be charged under new, currently unknown gauge transformations. Such 
interactions, if at all present, are neglected in our analysis. 

The most general renormalizable Lagrangian consistent with SM gauge invariance is 

Cu = £oid - X m L a HN l - ~yN"n 1 + H.c, (II.l) 

z=l 

where C id is the traditional SM Lagrangian, H is the Higgs weak doublet, L a , a = e, ji, t, are lepton weak doublets, 
\ ai are neutrino Yukawa couplings, and Mj are Majorana masses for the N^, We choose, without loss of generality, 
the Majorana mass matrix Mr to be diagonal and its eigenvalues Mi to be real and positive. We also choose the 
charged lepton Yukawa interactions and the charged weak current interactions diagonal so that all physical mixing 
elements are contained in the neutrino sector. 

After electroweak symmetry breaking (when H develops a vacuum expectation value v), C u will describe, aside 
from all other SM degrees of freedom, six neutral massive Weyl fermions — six neutrinos. The resulting mass terms 
can be expressed as: 




(II.2) 



where /i = Xv, and v and N are vectors of the three active neutrinos {v e , v^, v T ) and the three right-handed, sterile 
states, respectively. Each entry in the symmetric mass matrix of Eq (jll.2(l is itself a 3 x 3 matrix of mass parameters. 
Diagonalization of the mass matrix yields eigenstates with Majorana masses that mix the active-active states, related 
via the standard lepton mixing matrix V , and the active-sterile states. The physical neutrinos are thus linear 
combinations of all active and sterile states. Throughout, we will work in the "seesaw limit," defined by \i <C Mr. 
In this case, there are three mostly active light neutrinos and three mostly sterile heavy neutrinos where 'mostly' is 
determined by the induced mixing parameters. 

In the seesaw limit, the diagonalization is simple, are real: 

( ^ k ) - ( -\A ) ( l ! ) ( o l ) ( ^ l ) ( ■? ) + °(^. <n-3) 

where m is the diagonal matrix of light neutrino masses and is a matrix of active-sterile mixing angles found from 
the relations 

QMr = fx, (II.4) 
QM R e T = -VmV T . (II.5) 

The elements of Q are small (0(h/Mr)), and the standard seesaw relation (VmV T = — fiM^ fi T ) is easily obtained 
by combining Eqs. (|II.4(I and Eqs. I|II.5(1 . On the other hand, using Eq. I|H.5|1 . we can relate the mixing parameters in 
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to the active-active mixing angles contained in V and the neutrino mass eigenvalues. In the case of observably light 
sterile neutrino masses, as considered in our analysis, this equation is very useful, as it places testable constraints on 
observable quantities. The general solution (first discussed in detail in Q) of Eq. 1)11.5(1 is 

9 = -Vm 1/2 OM R 1/2 , (11.6) 

where O is an orthogonal matrix parameterized by three mixing angles 0i2, </>i3, <?^23-* Physically, the mixing matrix 

is a consequence of our freedom to choose the form of Mr. An illustrative example is found by considering the 
mass ordering of Mr in its diagonal form. The reordered matrix M^(mj <-> rrij) is equivalent to a 7r/2 rotation in the 

1 — j plane and therefore represents the same physics as the original matrix, as it should. In other words, OMrO t is 
the physically relevant object, as opposed to O and Mr separately. This object, when constrained to be real, has six 
free parameters that we shall refer to as "heavy parameters:" (f)\2, <j>i3, (f>23, Mi, M2, M3. 

Using Eq. 1(11. 6J) . the 6x6 unitary neutrino mixing matrix is 

Note that, up to 0(0 2 ) corrections, V is unitary. U is entirely described in terms of the three light mixing angles, 
six mass eigenvalues and three angles <pij- Many combinations of these have been measured or constrained via 
oscillation searches, cosmology and astrophysics. In particular, the two active neutrino mass-squared differences and 
mixing angles have been measured Q, 0, [H| , thus leaving free the six heavy parameters and the absolute scale of 
active neutrino masses.' 1 ' With U, the corresponding neutrino mass values and the SM couplings we can calculate all 
observable quantities and compare them with data. 

It is natural to wonder how well the seesaw approximation holds once one starts to deal with Mr values around 
1 eV. From Eq. (|II.4|I , it is clear that one can choose for the expansion parameter yj tii/Mr. In all scenarios considered 
here, \JrajMR < 0.5 (for M ~ 1 eV and m ~ 0.3 eV). In the worst case scenario, therefore, first order corrections 
are 55% of the leading order terms, while second order corrections are near 30%. Corrections to most observables 
of interest are much smaller than this because they are suppressed by larger right-handed neutrino masses. The 
first non-trivial correction to Eq. (|II.6|) occurs at second order and we find that, for the ambitions of this paper, all 
approximations are under control. This simple argument has been verified numerically for the most worrisome cases. 

Before concluding this section, we wish to add that operators that lead, after electroweak symmetry breaking, to 
Majorana masses for the left-handed neutrino states (Mr,) are also allowed if one introduces SU(2)l Higgs boson 
triplets or nonrenormalizable operators to the SM Lagrangian. While we neglect these "active" Majorana masses, we 
caution the reader that the existence of such terms would alter our results significantly. In particular, assuming the 
seesaw approximation holds, Eq. I|II.5|) would read 

VmV T + OM R T = M L , (II.8) 

which leads to a relationship between 0, m, and Mr different from Eq. (|II.6|I . If this were the case, for ex- 
ample, it would no longer be true that the largest value (in absolute value) is constrained to be smaller than 
(^maxM/R.mm) 1 / 2 , where m max is the largest element of m, while Airunin is the smallest element of Mr. On the 
other hand, all objects on the left-hand side of Eq. I|II.8|) are observables. Hence, in the case of a low-energy seesaw, 
one can expect, in principle, to be able to test whether there are contributions to the neutrino mass matrix that are 
unrelated to the presence of right-handed neutrinos. By measuring V, m, Mr, and 0, one can establish whether Ml 
is consistent with zero. 



III. OSCILLATION PHENOMENOLOGY AND CURRENT EVIDENCE FOR LOW-ENERGY SEESAW 

Here we examine a number of experimental and observational anomalies that may be explained by light sterile 
neutrinos. More specifically, we explore what these can teach us about the currently unknown parameters of the 
seesaw Lagrangian, described in detail in Sec.^ In all cases we assume 3 mostly active and 3 mostly sterile neutrinos 
and, most of the time, will concentrate on a 3 + 2 + 1 picture of neutrino mass eigenstates, that is, three mostly active 



* In general, O is a complex orthogonal matrix. Here, however, we will restrict our analysis to real neutrino mass matrices, unless otherwise 
noted. 

1 The active neutrino mass hierarchy, normal vs. inverted, is another (discrete) free parameter. 
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sub-eV neutrinos, two mostly sterile eV neutrinos and one almost completely sterile keV neutrino. The hope is that 
the heavier state can account for warm dark matter (section llll Bll or pulsar kicks ( Sec. 1111 which both require at 
least one keV neutrino, while the other two mostly sterile states help "explain" the existing oscillation data where, 
for all practical purposes, the heaviest neutrino decouples and we are left with an effective 3 + 1 or 3 + 2 picture. We 
remind readers that a third possibility (2+2*) is currently ruled out by solar and atmospheric data [l2l \13L Il4| and 
will be ignored. 3 + 1 schemes that address the LSND anomaly are also disfavored by global analysis of short baseline 
oscillation experiments 0, 0, 0, 0] , and for this reason, we mostly concentrate on 3 + 2 fits to the LSND anomaly 

Our analysis method is as follows: For each experimental probe considered we perform a x 2 "fit" of the mixing 
matrix U, given by Eq. H11.7fl . and neutrino masses to the "data", and extract the region of parameter space that 
best explains the data. In most cases, we allow the light mixing angles and mass squared differences to vary within 
their ler limits (according to [To|]').§ the angles (f>i to vary unconstrained within their physical limits of — 2tt, and the 
lightest active neutrino mass eigenvalue to; to vary unconstrained between — 0.5 eV. The quotation marks around 
"fit" and "data" are meant to indicate that our methods are crude, in the sense that we are fitting to previously 
processed experimental data, assuming a diagonal correlation matrix with Gaussian uncertainties. In order to avoid 
the subtleties involved in such a "fit to a fit" , we hesitate to mention actual confidence intervals, but are compelled 
to do so for lack of a better measure of an allowed region. We sometimes present our best fit parameter points along 
with confidence intervals, but warn the reader to avoid strict interpretations of these numbers. While crude, our 
methodology of error analysis and fitting provides a very useful instrument for identifying whether (and how) the 
low-energy seesaw can accommodate a particular combination of data sets. 

Before proceeding, it is useful to cement our notation. The neutrino masses will be ordered in ascending order of 
magnitude from mi to me in the case of a normal active neutrino mass hierarchy (m\—m\ < m\—m\), while in the case 
of an inverted mass hierarchy they are ordered m% < mi < m2 < m^ < m§ < m§ (in this case |m| — m\\ > m\ — m\ ). 
The states with masses mi,2,3 are mostly active, while those with masses 7714,5,6 are mostly sterile. Elements of the 
mixing matrix are referred to as U a %, where a = e, fx, t, si,S2, S3 (s's are the right-handed neutrino degrees of freedom) 
and i = 1,2, 3, 4, 5, 6. We also define Am^ = mj — to 2 and will refer to the lightest active neutrino mass as to;. In 
the case of normal (inverted) active neutrino mass hierarchy m; = toi (to; = 7713). 

A. Short baseline oscillation constraints 

Here we analyze the constraints imposed on the unknown mixing parameters by current neutrino oscillation data. We 
will assume that all solar, reactor, long-baseline and atmospheric data are properly fit with active-active oscillations, 
and that constraints on the other seesaw parameters will be provided mostly by short-baseline accelerator experiments. 
It is interesting to note that the inclusion of the angles introduces enormous freedom into the system. Any one 
active-sterile mixing angle contained in O can always be set to zero by an appropriate choice of O. In fact, all but 
three elements may be set to zero simultaneously, with only a single non-zero element in each row and column. In this 
case, these are constrained to be around y to/ /Mm where i is the column of the non-zero element. This is especially 
true when the mostly active neutrino masses are quasi-degenerate. An important "sum rule of thumb" is the following. 
For a given right-handed neutrino mass Mj, the active-sterile mixing angle squared is of order m/Mi, where to is a 
typical active neutrino mass. One can always choose parameters so that, for at most, two values of a = e,ji,T, U a i 
are abnormally small. In that case, however, the "other" U a i is constrained to saturate the bound |C/ Q i| 2 < mi/Mi. 

The most compelling evidence for light sterile neutrinos comes from the short baseline oscillation experiment by 
the Liquid Scintillator Neutrino Detector (LSND) collaboration at Los Alamos. Using a ~ 30 MeV beam they 
observed a better than 3cr excess of i7 e -like events above their expected background at their detector some 30 m 
away from the production point 0. This evidence of V e oscillation requires a mass-squared difference greater 

than 1 eV 2 , clearly incompatible with the small mass-squared differences observed between the active neutrinos. 
Several mechanisms, such as CPT- violation 0, Lorentz invariance violation ^t|> quantum decoherence [l8| . 
sterile neutrino decay |l9j| and, of course, oscillation into sterile neutrinos have been proposed to explain this result. 
Here we concentrate on the last possibility. 

In order to take into account all short baseline data we "fit" our mixing parameters and masses to the results of 
the 3 + 2 performed in 45], which are summarized in Table ITI |20j| . Here we assume that the heaviest, mostly sterile 
state does not participate effectively in LSND oscillations. This is guaranteed to happen if me > 10 eV. On the other 



t It would have been rather difficult to construct a 2+2 neutrino mass hierarchy using the seesaw Lagrangian. 

§ In the case of 3+1 "fits" (c/. Eqs. 1111.101111.111 1. we kept the active neutrino parameters fixed at their best-fit values. 
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hand, |J7 Q 6| 2 are partially constrained by our attempts to accommodate LSND data with seesaw sterile neutrinos, as 
will become clear in the next subsections. 

We find that mi, the lightest neutrino mass, is constrained to lie between (0.22 — 0.37) eV, with a "best fit" value 
of 0.29 eV. Thus, the active neutrino mass spectrum is predicted to be quasi-degenerate. A sample 6x6 neutrino 
mixing matrix that fits all data is 



U3+2 — 



/ 0.8301 


0.5571 


0.001365 


0.1193 


-0.009399 


-0.006513 \ 


-0.3946 


0.5866 


0.7072 


0.2016 


0.2262 


0.0003363 


0.3932 


-0.5879 


0.7070 


0.4760 


-0.0949 


0.001470 


-0.2067 


0.09514 


-0.4792 


1 








0.1343 


-0.1832 


-0.09284 





1 





\ 0.004963 


0.004295 


-0.001268 








1 / 



(HI.9) 



while the associated masses are mi ~ 7772 — «^3 = 0.28 eV, 7774 = 1.0 eV, 777,5 = 4.7 eV, and me — 6.4 keV. Note that 
the matrix in Eq. (|III.9ll is only approximately unitary, up to corrections of order 25%. This result agrees qualitative 
with those obtained in [3j. The neutrino masses and mixings obtained in this "fit" are depicted in Fig.Q] We will use 
the results of "fits" similar to this one throughout the paper. 

One can also aim at a (currently disfavored) 3 + 1 LSND fit.^. In this case, much lower mi values are also allowed, 
extending well into the hierarchical spectrum range. In this case, all mi values above 0.01 eV and 0.03 eV are allowed, 
assuming an inverted and normal mass hierarchy, respectively. This is to be compared with the results found in 0], 
where only trivial choices for O were considered. Examples that "fit" all oscillation data include, for an inverted 
active mass hierarchy: 7774 ~ m-2 — 0.066 eV, 7773 = 0.043 eV, 7774 = 0.96 eV, 777,5 = 5 keV, and 777g = 10 GeV, together 
with 



U. 



rtcd 



3+1 



0.8305 
-0.3939 

0.3939 
-0.01486 
0.001370 



0.5571 
0.5872 
-0.5872 
-0.2218 
-0.001878 



\ 0.000002372 0.0000004094 





0.7071 
0.7071 
-0.1134 
0.002253 
-0.0000007187 



0.1359 
0.2046 
-0.04421 

1 







-0.00009142 
0.00005000 
-0.003236 



1 





-0.000002198 \ 
0.000001202 
-0.0000001857 





1 / 



(111.10) 



For a normal mass hierarchy, we find that mi 
7775 = 5 keV, mg = 10 GeV, and 



0.055 eV, 777 2 = 0.056 eV , 7773 = 0.0744 eV, m 4 = 0.96 eV, 



U. 



normal 
3+1 



0.8305 
-0.3939 

0.3939 
-0.05028 
0.0008214 



0.5571 
0.5872 
-0.5872 
-0.1355 
0.002545 





0.7071 
0.7071 
-0.2231 
-0.002310 



\ -0.000002220 0.0000007646 0.00000005663 



0.1173 
0.2176 
0.09802 

1 







-0.002100 
0.0004625 
0.002804 



1 





0.000001418 \ 
-0.000001364 
0.000001283 





1 / 



(111.11) 



"fit" all oscillation data quite well. 

Note that a null result from MiniBooNE is bound to place significant limits on the seesaw energy scale. If all right 
handed neutrino masses are similar, the effective mixing angle that governs za, — » v e transitions is sin 2 26*MiniBooNE 



< 



TABLE I: Parameter values used in our analysis. These were extracted from a fit to all short baseline neutrino oscillation 
experiments including LSND within the 3 + 2 scenario [Tnl l2(l | . la indicates a rough estimate of the 1 sigma allowed range for 
the different parameters. 





u e4 




u e5 


Uei_i5 


Amli (eV 2 ) 


AU& (eV 2 ) 


Central Value 


0.121 


0.204 


0.036 


0.224 


0.92 


22 


la 


0.015 


0.027 


0.034 


0.018 


0.08 


2.4 



■ This is easily accomplished by requiring 7715 > 10 eV. 
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FIG. 1: Neutrino mass eigenstate spectrum, along with the flavor composition of each state. This case accommodates all 
neutrino oscillation data, constraints from r-process nucleosynthesis in supernovae, and may help explain anomalous pulsar 
kicks (see text for details). While we choose to depict a normal hierarchy for the active neutrino states, an inverted active 
neutrino mass hierarchy would have yielded exactly the same physics (as far as the observables considered are concerned). 



4m 2 /M 2 . Hence, a null result at LSND would rule out a seesaw energy scale M lighter than 6 eV, assuming all 
active neutrino masses m are around 0.1 eV Q- This limit is sensitive to the lightest neutrino mass mi and can be 
somewhat relaxed (similar to how we obtain a good 3+2 to all neutrino data) by postulating a (mild) hierarchy of 
right-handed neutrino masses and by assuming that sterile-electron and sterile-muon neutrino mixing is suppressed 
with respect to naive expectations for the lightest mostly sterile state(s). For larger values of mi, M values around 
10 eV are already constrained by — > v e searches at the NuTeV j2l| and NOMAD 22] experiments, and v„ — > v T 
searches at CHORUS 0. 



B. Cosmological and Astrophysical Constraints, Warm Dark Matter 

Very light sterile neutrinos that mix with the active neutrinos are constrained by several cosmological and astro- 
physical observables. The "seesaw" right-handed neutrinos are no exception. Given that active-sterile mixing angles 
| C/o-i | 2 < mi/mi (a = e, p, r, i = 4, 5, 6), it turns out that for "standard cosmology," the right-handed neutrinos ther- 
malize with the early universe thermal bath of SM particles, as long as the reheat temperature is higher than their 
Majorana masses. For the low seesaw energy scales we are interested in, this is a problem. For the values of Mr under 
consideration here, thermal right-handed neutrinos easily overdose the universe. Smaller mi values [mi < 10 -5 eV) 
lead to the possibility that right-handed neutrinos are the dark matter, as recently discussed in the literature 0,0- 
We comment on this and other possibilities shortly. 

Fig. [21 depicts the region of the |?7 e 6| 2 x mg-plane in which the contribution of the heaviest neutrino vq to Q, (the 
normalized energy density of the universe, p/p c ) is larger than 0.3 (dark region). The same constraint roughly applies 
for all a = e, p,r and i = 4, 5, 6. The dashed diagonal lines correspond to |C/ Q g| 2 = mi/me, for different values of mi. 
All lines lie deep within the dark Sl s > 0.3 region. 

For smaller values of Mr, the situation is also constrained. For Mr values below tens of eV, thermal sterile neutrinos 
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FIG. 2: Adapted from |4fll |. Cosmological and astrophysical constraints on the j[/ Q 6| 2 X m6-plane. In the large dark grey region, 
the density of a thermal v§ population is f2 s > 0.3, while the light grey 'X-ray' region is disfavored by X-ray observations. The 
regions labeled 1,2,3 are preferred if one is to explain anomalous pulsar kicks with active-sterile oscillations inside supernovae. 
Regions 1 and 3 qualitatively extend inside the f2 s > 0.3 part of the plane as indicated by the horizontal dotted and dash- 
dotted lines, respectively. The regions 'Warm Dark Matter" and "Too Warm Dark Matter" are meant to represent the region 
of parameter space where thermal v§ qualifies as a good (or bad) warm dark matter candidate. The region above the solid 
diagonal line is disfavored by the observation of electron (anti)neutrinos from SN1987A. The diagonal dashed lines correspond 
to !7e6 = mi /ma, for different values of mi. Also shown is our "best fit" sterile solution for different pulsar kick scenarios, 
assuming the 3 + 2 LSND fit for the lighter states. The regions one and three best fit values are represented by circles and a 
star, respectively. See text for details. 



contribute to the amount of hot dark matter in the universe |2J, l25l 1261 ] . Right-handed neutrinos will thermalize as 
long as m, sin 2 8 la > 5 x 10" 4 eV H H|. In low energy seesaws, this roughly translates into to > 10 3 eV, where 
to is the active neutrino mass scale. For the cases of interest here, to > yj Ato 2 ^ ~ 10~ 2 eV, in which case to^ values 
above 0.6 eV are ruled out (according to 26]).* Note that any sterile neutrino solution to the LSND anomaly faces a 
similar problem, which must be resolved with non-standard cosmology. 

One way to avoid the bounds described above (see also, for example, [s0|) is to consider that the reheating tem- 
perature of the universe is very low, say T r < 5 MeV |3l|. This way, right-handed neutrinos, in spite of their "large" 
mixing angles, never reach thermal equilibrium in the early universe and neither overdose the universe, nor contribute 
to the amount of hot dark matter. Unless otherwise noted, this is the assumption we make here. 

Under these circumstances, it is interesting to consider whether light seesaw right-handed neutrinos still qualify as 
good warm dark matter. This could happen if their production in the early universe was non-thermal. One concern 
surrounding warm dark matter is whether it is ruled out by large scale structure surveys. Here, we will not add to 
this discussion but refer readers to the recent literature on the subject j33|. A brief summary of the situation is as 
follows: constraints on warm dark matter can be translated into a lower bound on the mostly sterile neutrino mass. 



* In our "best fit" 3+2 solution to the LSND data, the sum of all active neutrino masses violates slightly the constraint obtained in |2(I . 
This problem can be easily evaded if we choose mj values close to the lower end of the "allowed region" . Other "mild" non-standard 
cosmology effects f2g| are also known to alleviate the hot dark matter bound on neutrino masses. 
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The lower bound has been computed by different groups, and lies somewhere between 3 and 14 keV 32] . Different 
lower bounds depend on several issues, including which subset of Lyman alpha-forest data is taken into account. 

Another constraint on potential dark matter sterile neutrinos comes from the observation of X-rays originating in 
galactic clusters. Such regions of the universe should be overdense with warm dark matter heavy neutrinos, which can 
be directly observed via their radiative decay v§ — > Vi +7 |33|. Bounds from X-ray observations have been summarized 
very recently in [^H • Combining the results of j^S] and Fig. |3 we find that for lightest neutrino masses larger than 
10~ 2 eV, such bounds can only be avoided for m, < 100 eV, where large scale structure constraints on warm/hot 
dark matter are severe. This qualitative analysis indicates that seesaw sterile neutrinos cannot simultaneously fit the 
LSND data and serve as cold dark matter. 

On the astrophysics side, the most severe constraint on light, sterile neutrinos is provided by the observation of 
electron (anti)neutrinos coming from SN1987A. The current analysis consists of comparing the model-dependent 
neutrino flux at the surface of the neutrinosphere with that detected on Earth. Large sterile neutrino mixin g an d 
mass would result in modification/depletion of the detected neutrino signal (for a recent detailed discussion, see 35]). 
Although only twenty neutrinos were observed in this event, one can still place bounds on sterile-active neutrino 
mixing. As far as "LSND" sterile neutrinos are concerned, these bounds are still weaker than those obtained by the 
null short baseline oscillation experiments |36| and therefore already accounted for in our analysis. Heavier right- 
handed neutrinos can, however, be excluded by SN1987A neutrino data. Fig. |2] depicts the region of parameter space 
excluded by SN1987A data (region to the right of solid, diagonal line). This bound is defined by mi\pl\/X) r a i >0.22 keV 
[3ll l37| ] . See also [28| . According to Fig. |3 supernova bounds force the seesaw scale to be below a few keV for to,/ 
values above 0.01 eV. 



C. Pulsar Kicks 

Pulsars are born from the gravitational collapse of the iron core of a massive star. These core collapse supernova 
are an excellent source of neutrinos, producing all (active) flavors copiously (see [3^ for a detailed review). Current 
observations point to the fact that some pulsars move with peculiar velocities much greater than those expected from 
an asymmetric supernova explosion mechanism. Quantitatively, current three dimensional models yield velocities up 
to 200 km/s j3s|, while pulsars moving at speeds as high as 1600 km/s have been observed. 

Since roughly 99% of the approximately 10 53 ergs of energy released in a core collapse supernova is in the form of 
neutrinos, it is reasonable that neutrino physics provides a solution to this anomaly. At these rates a small (1 — 3)% 
asymmetry in neutrino emission can account for the observed large pulsar velocities. Neutrinos are always produced 
asymmetrically in the polarized medium of the proto-neutron star, due to the left-handed nature of their interactions. 
Unfortunately, asymmetric production cannot solve this problem because the associated medium densities are such 
that neutrinos undergo multiple scattering within the star's interior, eventually diffusing out of an effective surface, 
called the neutrinosphere, with all initial asymmetries washed away. Several distinct mechanisms have been formulated 
to sidestep this fact. Specifically, the existence of large neutrino magnetic moments has been explored in [4(i|. and 
can be tested in next generation neutrino scattering experiments ^], 0; E3 • Proposed solutions also exist which 
utilize standard three flavor neutrino oscillations, where and z/ T , appearing between their neutrinosphere and the 
larger v e neutrinosphere, can stream unhindered out of the star |44|. This solution is, however, currently disfavored 
by terrestrial oscillation experiments. 

We concentrate on the case of oscillations into sterile neutrinos, which can proceed in various ways, depending on 
the mass and coupling of the relevant neutrinos as well as the properties of the collapsing star, including its density 
and magnetic field. Following [4^, we separate and analyze these within three distinct categories. Each one requires 
the existence of a keV-scale sterile neutrino with very small couplings to the active flavors, of the order 10~ 4 — 10 -5 , 
especially if light sterile neutrinos are thermally produced in the early universe. Under these circumstances, if seesaw 
neutrinos are to play the role of the sterile neutrinos responsible for pulsar kicks, |£/Q,j| 2 < m/rrii [i = 4, 5, 6) must 
lie in the 10 -9 range for m t ~ 10 4 eV. This implies to ~ 10 -5 eV and is only compatible with a hierarchical active 
neutrino mass spectrum and very light to,/ , as identified in 0, @ • 

Here, instead, we will concentrate on identifying solutions that will address pulsar kicks and the LSND anomaly. 
According to the discussion in the previous subsection, the mostly sterile neutrino masses TO4 and TO5 are con- 
strained to be less than 10 eV so that a 3+2 solution to the LSND anomaly can be obtained from the seesaw 
Lagrangian. The heaviest neutrino mass mg is unconstrained, so we are free to vary it as needed in order to attack 
the pulsar peculiar velocity issueJ Naively, the fraction of v a (a — e,/i,r) in vq is expected to be of the order 



t We can neglect the lighter sterile neutrinos [1/4 and fg), as they should not alter the kicking mechanism significantly due to their small 
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U a 6 ~ v/0^3 eV/3 x 10 3 cV = 10~ 2 , much too large to satisfy the pulsar kick plus cosmology constraints summarized 
in Fig. [3 On the other hand, once the Jig < 0.3 constraint is removed, the 'pulsar kicks' allowed region of the plane is 
significantly enlarged, as qualitatively indicated by the horizontal lines in Fig. [21 In this case, which we must consider 
anyway if we are to have agreement between LSND and searches for hot dark matter, one can envision explaining 
pulsar kicks and the LSND data simultaneously. 

In scenario 1, the pulsar kick is produced via an active-sterile MSW resonance in the core of the proto-neutron 
star at large densities, greater than 10 14 g/cm 3 , and magnetic fields, near 10 16 G 0|. The effective neutrino matter 
potential in material polarized by a strong magnetic field contains a term proportional to k ■ B/\k\ |47L l48j . where k 
is the neutrino's three-momentum and B is the local magnetic field vector. Clearly the MSW resonance occurs at 
a radius that depends on k ■ B/\k\, the relative orientation of the neutrino momentum and magnetic field. Sterile 
neutrinos produced at smaller radii (higher temperatures) carry greater average momentum than those produced 
at larger radii (lower temperatures), yielding an asymmetric momentum distribution of emitted neutrinos. This 
asymmetry is capable of producing the observed pulsar kicks, in the direction of the magnetic field, when the mass 
and coupling of the sterile state is near 8 keV and above 1.5 x 10 -5 , respectively |4^. We found the "best fit" to 
the LSND data (using 1/4 and v^) and pulsar kicks (using v§ ) and mg > 5 keV. The |£/ a 6| 2 an( i m 6 "best fit" values 
are depicted in Fig. This solution is strongly disfavored by the observation of neutrinos from SN1987A. The fact 
that |J7 e g| is much larger than the other two active-sterile mixing angles is due to the fact that |J7 e 4| and \U e s\ are 
constrained by LSND data to be much smaller than naive expectations (see Eq. (|III.9|) '). In order to reduce |U e 6|, one 
would have to either reduce mi by an order of magnitude - which renders the 3+2 fit to oscillation data very poor - 
or increase me, which would only push \U e e\ deeper into the region of parameter space ruled out by SN1987A. One 
can, however, find 3+1 solutions to LSND data where could pose as the sterile neutrino that explains why pulsar 
peculiar velocities are so large (see Eqs. (|III.10IIII.lip L 

Scenario 2 also relies on a direction-dependent MSW resonance, this time occurring outside the core where the 
matter density and temperature are much lower. Here, both the active and sterile neutrinos are free to stream out 
of the star. The departing active flavors still have a small interaction cross-section, a ~ G^-E 2 , and can therefore 
deposit energy and momentum into the star's gravitationally bound envelope proportional to the matter it transverses. 
Via the direction-dependent resonance, neutrinos moving in the direction of the magnetic field remain active longer, 
deposit more momentum, and thus, kick the star forward. The observations can be explained in this case with a 
smaller sterile neutrino mass and larger active-sterile coupling near 4 keV and 4.5 x 10~ 5 , respectively |4^. In the 
case of our LSND "fit" the data, we can constrain one of 1 2 or |£/ T 6| 2 to lie inside region 2. The other |J7 Q g| 2 
(a = e, t or e, /1), however, are constrained to be large, and lie inside region 3. 

Scenario 3 proceeds through off-resonance production of the sterile neutrino in the proto-neutron star core |49| . 
The amplitude for sterile neutrino production by a weak process is proportional to U£q, the effective mixing angle 
between the heavy mass eigenstate and the flavor eigenstate. Initially, this quantity is very small due to matter effects 
in the dense core. The effective potential in the star's interior is quickly driven to zero in the presence of sterile 
neutrino production by a negative feedback mechanism. If this occurs in a time less than the diffusion time-scale for 
the active neutrinos, approximately (3— 10) s, the mixing angle will reach its vacuum value |50|. The sterile neutrinos 
will then be produced and emitted asymmetrically and thus, kick the pulsar to large velocities. Lower limits on the 
vacuum mass and mixing values are derived by requiring that the off-resonance time scale (inversely proportional to 
m| sin 2 26 a o) for the evolution of the matter potential to zero be less than about ten seconds. This places the sterile 
mass and mixing at approximately 1 keV and above 5 x 10~ 5 , respectively. Since all three active flavors are present 
in equal abundances and all contribute to the effective matter potential, the mixing angle in question is not any 
particular U a Q. Rather, it is the angle 9q associated with the projection of v§ onto the space spanned by v ei and 
u T) that is Q\ = [/ 2 6 + £/ 2 6 + U^ 6 . From Eq. (|II.6(1 we see that 9§ = \fmij m§ up to corrections due to the non-unitarity 
of V and active neutrino mass differences. This is independent of mixing angles, and therefore cannot be tuned to be 
small. Our "best fit" region-3 solution is depicted in Fig. by a star. It turns out that U a § have very similar values 
for a = e,fi,T. In order to evade the SN1987A constraint, we were forced to pick m; values close to lower bound of 
our 3+2 LSND "fit" (mi = 0.22 eV), so that \U e 4\, {U^], and | 1 are close to the low-end of the allowed range in 
TableHI 



mass and non-resonant production. 
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D. Supernova Nucleosynthesis 

Core collapse supernova are believed to produce the observed heavy element (A > 100) abundance through the r- 
process, or rapid neutron capture process. Here we briefly review this mechanism (see~|5l| for a comprehensive review), 
as well as its facilitation by the addition of active-sterile neutrino oscillations |§2 ■ This scenario begins in the neutrino 
driven wind; that is, the wind of ejected nucleons driven by neutrinos radiated from the cooling proto-neutron star. 
The maintenance of equilibrium among neutrons, protons, and electron (anti)neutrinos in neutrino capture processes 
leads to a neutrino-rich environment. As the wind propagates, it cools enough for all free protons to bind into a 
particles. In the ideal r-process picture, as the wind cools further these a particles bind into intermediate size seed 
nuclei which later undergo neutron capture to form the observed heavy r-process elements. 

This ideal scenario is dampened by the large number of electron neutrinos present at the stage of a particle 
formation. These will capture on the free neutrons, converting them to protons, which in turn will fuse to make 
more a particles. The end result is a very small free neutron to a particle ratio, conditions unfavorable for r-process 
element formation. This is known as the a effect and must be circumvented to produce the correct distribution of 
heavy elements. A clear solution to this problem is to reduce the number of electron neutrinos present at this stage, 
which can be accomplished by resonant v e — + v s conversion^ |52| . 

The sterile neutrino solution to the r-process mechanism is modeled and fit to the data in reference [E^ including the 
effects of relevant nuclear physics and additional neutrino oscillations in the star's envelope. The analysis is expanded 
in 01 with the inclusion of fission cycling of the produced heavy elements. The analysis indicates the need for an 
eV-scale sterile neutrino with an allowed parameter space much larger than that constrained by LSND. By itself, the 
requirement of successful r-process in supernovae only weakly constrains the light neutrino mass scale to be greater 
than 10~ 2 eV and 10 -3 eV for a 1 eV and 10 eV sterile neutrino, respectively. With regard to the LSND results, 
it has been demonstrated that the 3+1 oscillation scenario fits within this parameter space [54j . Considering that 
the best fit mass-squared difference and mixing angles for the fourth mass eigenstate, which makes up most of the 
lightest sterile neutrino, is very similar between the 3 + 1 and 3 + 2 case |20|. it is reasonable to conclude that v e <-» v S i 
resonant conversion will also fit within this scenario. Even oscillations into the heavier is s5 state can potentially solve 
this anomaly if the neutrino driven wind expansion time-scale is sufficiently small, < 0.1 sec. To conclude this section 
we note that, although the sterile neutrino solution to the supernova nucleosynthesis problem fits well within our 
seesaw framework, it adds no additional constraints, and therefore does not increase the predictability of our scenario. 

IV. OTHER PROBES OF THE SEESAW ENERGY SCALE 

Here we survey other existing and future probes of light sterile neutrinos. As opposed to the previous cases, these 
probes are perfectly consistent by themselves. That is, extra heavy neutrinos are not required to solve problems 
within the system. However, their addition can lead to large modifications to the outcome of such experiments, 
thus rendering the eV-scale seesaw scenario testable. Specifically, we consider bounds from tritium beta-decay and 
neutrinoless double-beta decay. Observations in all of these areas have already yielded useful constraints on sterile 
neutrinos, and the situation is expected to improve in the next few years. 

A. Tritium Beta-Decay 

The endpoint of the electron energy spectrum in the beta-decay of tritium is a powerful probe of nonzero neutrino 
masses. This results from the decay kinematics of the system, which is necessarily modified by the presence of a 
massive neutrino. The nonzero neutrino mass effect can be understood almost entirely from the analysis of the phase 
space distribution of the emitted electrons, and is therefore quite model independent. Existing beta decay experiments 
extract limits on an effective electron neutrino mass m 2 = \ U e i\ 2 mf [55j, provided that the neutrino masses are 
smaller than the detector energy resolution. Currently the most stringent bounds on m 2 are (2.3) 2 eV 2 at 95% 
confidence from the Mainz experiment [5(j and (2.5) 2 eV 2 at 95% confidence from the Troitsk experiment [H^. In the 
next few years the Katrin experiment should exceed these limits by nearly two orders of magnitude, probing down to 
(0.2) 2 eV 2 at the 90% confidence level [58|. One might naively compute this effective mass for the "best fit" mixing 



+ In this mechanism the effective matter potential, which depends on the number of electrons, positrons and neutrinos, varies wildly as a 
function of distance from the core. Along this radial dire ction there are three relevant MSW resonant conversions that must be tracked 
and understood: u e — * u s , v e — > u s and v s — » v e . See |52| for more information. 
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FIG. 3: 1 — S/So as a function of the /3-ray energy, where S is the /3-ray energy spectrum obtained assuming three mostly 
active, degenerate neutrinos with mass m — 0.1 eV and one mostly sterile neutrino Ui with mi = 0.1, 0.5, 1, 5, and 10 eV. The 
mixing angle is given by = m/rrii. So is the spectrum associated with massless neutrinos. See text for details. 



parameters obtained in the previous section. In this case, we expect a keV seesaw neutrino to contribute to m 2 by 
a huge amount dm';, = U^m 2 - - - ■•'• 2 



'e6'"6 m6 
-3 



— m\ = mime [3| so that it would be excluded by the current precision measures 



of tritium beta-decay for mi > 10 eV. This is clearly an incorrect treatment of the physics. As pointed out in, for 
example, 59], the existence of a heavy neutrino state would produce a kink in the electron energy spectrum of size 
\U e i\ 2 at an energy Eo — mi as well as a suppression of events at the endpoint of order 1 — |t/ e i| 2 . Here f7 e , (i = 4, 5, 6) 
is the mixing between the electron neutrino flavor eigenstate and the heavy mass eigenstate, while Eq — 18.6 keV is 
the endpoint energy of tritium beta-decay. 

Fig. 01 depicts 1 — S/Sq, where S is the /3-ray energy spectrum obtained assuming three mostly active, degenerate 
neutrinos with mass m = 0.1 eV and one mostly sterile neutrino Ui with various masses m, and mixing angle 
— m/rrii, while So is the spectrum associated with massless neutrinos. One can readily observe "kinks" in the 
spectrum above m,. For /3-energies above Eq — mi, the impact of the sterile state is to "remove" around 1 — \U e i\ 2 of 
the /3-rays from spectrum. This is most significant between Eq — mi and Eq minus the mass of the active neutrinos. 
For energies below Eq — mi, the spectrum agrees with that obtained from the emission of one effective neutrino with 
mass-squared m 2 c . 

We estimate the sensitivity of future tritium beta-decay experiments to the emission of one heavy state by con- 
sidering the ratio between the number of electrons with energies above Eq — AE in the case of one heavy massive 
neutrino Vi and in the case of massless neutrinos, 



R(U ei ,M R ) 



\U K 



E -AE 



dE^L( mi ) + (1 - \U el \ 2 ) 



Eq-AE 




En-AE 



<<»> • 



(IV.12) 



where dN/dE is the energy distribution of /3-rays, which depends on the neutrino mass rrii- This expression can 
be easily generalized for more than one heavy neutrino. The advantage of using the ratio above is that potential 
systematic uncertainties and normalization effects can be safely ignored. An experiment is sensitive to a massive 
neutrino state if it can distinguish R from unity, a determination that should be limited by statistics due to the very 
low /3-ray flux in the high-energy tail of the electron spectrum. 
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FIG. 4: Contour plot of constant R, as defined by Eq. 1IV.12L assuming an energy window AE = 25 eV. The solid (dashed) 
line corresponds to \J mi/mi, a naive upper bound for \U e i\, for mi — 0.3 eV (0.01 eV). The circles correspond to U e i for the 
three mostly sterile states obtained by our "fit" to other neutrino data, Eq. 1)111.90 . See text for details. 



In order to compute R, we use an analytic expression for Eq. I|1V.12|) . which exists provided that one neglects 
nucleon recoil in the decay. Fig. 01 depicts constant i?-contours in the \U e i\ x rrii plane. Contours were computed 
for AE — 25 eV, in order to allow one to easily compare our results with the sensitivity estimates of the Katrin 
experiment. After data-taking, Katrin is expected to measure R at the 0.1%-1% level (lightest grey region). Its 
sensitivity is expected to be ^/m^ > 0.2 eV at the 90% confidence level. This can be extracted from the plot by 
concentrating on the U e i — 1 line. Note that while the expected energy resolution for Katrin is of order 1 eV, the 
expected number of signal events above Eq — 1 eV is minuscule (both in absolute terms and compared with expected 
number of background events), so that most of the sensitivity to nonzero neutrino masses comes from analyzing the 
shape of the electron spectrum in the last tens of electron-volts. A larger "window" would suffer from increased 
systematic uncertainties, so that AE ~ 25 eV is representative of Katrin's optimal reach |58| . 

The shape of the constant R contours is easy to understand. As already discussed, for rrii > AE, the effect of the 
right-handed neutrinos is to reduce the spectrum in an energy independent way by 1— |J7 e i| 2 , while for rrii < states 
with the same effective mass-squared mf \U e i\ 2 produce the same effect in tritium beta-decay so that the diagonal 
lines coincide with lines of constant m 2 \U e i\ 2 . 

A more sensitive approach would be to "bin" the last tens of eV of the "data" into 1 eV bins, and fit the distribution 
to a massless neutrino hypothesis. For the values of the parameters in which we are interested, we find that one 25 eV 
bin yields roughly the same sensitivity to nonzero neutrino masses as twenty five 1 eV bins for large masses and small 
mixing angles. For smaller masses and larger mixing angles, a "binned" analysis should be sensitive to effects which 
are localized in individual bins (such as "kinks"). Another recent estimate of the sensitivity of Katrin to heavy, sterile 
neutrinos can be found in |28|. 

Fig. 0| also depicts the loose upper bound for U e i — 'mi /rrii as a function rrii, for mi = 0.32 eV and mi = 0.01 eV. 
For mi > 0.1 eV, Katrin should be sensitive to Mr < 1 keV, while for mi > 0.01 eV (the solar mass scale) Katrin 
should be sensitive to Mr > 10 eV and Mr < 100 eV, where here we assume that all right-handed neutrino masses are 
of order Mr. In the case of seesaw parameters that fit the LSND data with a 3+2 neutrino spectrum (see Eq. i|III.9|) ). 
expectations are high as far as observing a kinematical neutrino mass effect at Katrin, in spite of the fact that the 
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fit to LSND data requires U e 4 and f/ e 5 to be "abnormally" low. Fig. 0] depicts f7 e i and values for the heavy 
neutrinos (open circles). The contribution of the heaviest of the two LSND-related sterile neutrinos is of order the 
Katrin sensitivity, while the active contribution itself, which leads to m 2 e = J2i=i 2 3 \^ei m i \ ~ m f ^ s already within 
the Katrin sensitivity range, given that large mj > 0.22 eV values are required by our 3+2 LSND "fit". The effect of 
i>6 is small if me is larger than 1 keV (required if one takes the "pulsar kicks" hint into account), but would be very 
significant if ttiq were less than 1 keV. 



B. Neutrinoless Double-Beta Decay 



If the neutrinos are Majorana fermions - as predicted in the case of interest here - lepton number is no longer a 
conserved quantity. The best experimental probe of lepton number violation is the rate for neutrinoless double-beta 
decay. This process, which violates lepton number by two units, is currently the subject of intense search [(5(j. If 
neutrino masses are the only source of lepton number violation, the decay width for neutrinoless double-beta decay 
is 



-M(m 2 ,Q 2 ) 



(IV.13) 



where Ai is the relevant nuclear matrix element and Q 2 ~ 50 2 MeV 2 is the relevant momentum transfer. In the limit 
of very small neutrino masses (ra 2 <C Q 2 ), Fo u pp is proportional to an effective neutrino mass \m ee \, 

n 

m ee =J2 U !i m i- ( IV - 14 ) 

i 

The sum is over all light neutrino mass eigenstates. In the case of a low-energy seesaw, when all m t , i — 1, . . . , 6 are 
much smaller than Q 2 , it is easy to see that m ee vanishes The reason for this is that, in the weak basis we are 
working on (diagonal charged- lepton and charged weak-current), m ee is the ee-element of the neutrino mass matrix, 
as defined in Eq. (|II.2|) . One can trivially check that, by assumption, not only does m ee vanish, but so do all other 
m a /3, a, (3 = e, /1, r. Note that this result does not depend on the fact that we have been assuming all elements of the 
neutrino mass matrix to be real |6l| . 

For heavy z/j neutrinos, U^rrii no longer captures the dependency of To^pp on the exchange of z/j. For m 2 3> Q 2 , 
instead, the dependency on neutrino exchange is proportional to C/ 2 /to^. If this is the case, the overall contribution 
(including all heavy and light states) is no longer proportional to m ee but, instead, can be qualitatively expressed as 
a function of an effective m^f , 

The approximation T^pp cx \m^\ is good as long as one can neglect the dependency of M on m, and is not expected 
to be a great approximation when m 2 ~ Q 2 . Nonetheless, still qualitatively captures the behavior of To„pp as a 
function of the sterile neutrino masses and studying its behavior is sufficient for our ambitions in this discussion. 

Fig.[S]depicts for our "best fit" 3 + 2 LSND solution (see Sec lIII Ajl . as a function of the unconstrained me- As 
advertised, m|f vanishes for m\ -C Q 2 ■ The figure also depicts the "active only" value of m^ tlve = J2i=i 2 3 ^ei m i- 
Even in the limit m| 3> Q 2 , there is still partial cancelation between the mostly active and mostly sterile "LSND" 
states. This is a feature of the Lagrangian we are exploring here, and is not in general observed in other scenarios 
with light sterile neutrinos tailor-made to solve the LSND anomaly. 

Currently, the most stringent limits on this effective mass comes from the Heidelberg- Moscow experiment |62| where 
they find m ee < 0.91 eV at 99% confidence. In the near future, experiments aim to reach down to m ee values close to 
10~ 2 eV [6(j. A signal would rule out a seesaw scale below tens of MeV. On the other hand, if we were to conclude 
that the neutrino masses are quasi-degenerate (through, say, a signal in tritium beta-decay) and if the LSND 3+2 
solution were experimentally confirmed, a vanishing result for m ee could be considered strong evidence for a very 
small seesaw scale. On the other hand, if this were the case (m ee zero, large active neutrino masses), it would also 
be very reasonable to conclude that neutrinos are Dirac fermions. Distinguishing between the two possibilities would 
prove very challenging indeed. It is curious (but unfortunate) that in a low energy seesaw model, the neutrinos 
are Majorana fermions, but all "standard" lepton-number violating observables vanish, given that their rates are all 
effectively proportional to m a p\ 
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rn (eV): Heaviest sterile neutrino mass 

6 

FIG. 5: Effective m ee for neutrinoless double-beta decay as a function of me, the heaviest right-handed neutrino mass, assuming 
the existence of only light, active, neutrinos (magenta curve), with a degenerate mass spectra, and for our "best fit" 3 + 2 LSND 
sterile neutrino solution (blue curve). See text for details. Also indicated is the parameter region preferred by astrophysical 
hints of sterile neutrinos. We assume Q = 50 MeV. In the case of a low-energy seesaw, m ee vanishes as long as m§ <C Q. 



V. CONCLUDING REMARKS 

The "New Standard Model" (equal to the "old" Standard Model plus the addition of three gauge singlet Weyl 
fermions) is, arguably, the simplest extension of the SM capable of accommodating neutrino masses. This Lagrangian 
contains a new dimension-full parameter: the right-handed neutrino mass scale M, which must be determined ex- 
perimentally. Unlike the Higgs mass-squared parameter, all M values arc technically natural given that the global 
symmetry of the Lagrangian is enhanced in the limit M — > 0. 

Very large M values are theoretically very intriguing, and have received most of the attention of the particle physics 
community. There are several strong hints that new phenomena are expected at the electroweak breaking scale 
~ 10 3 GeV, the grand unified scale ~ 10 15-16 GeV, and the Planck scale ~ xq 18-19 GeV, and it is tempting to 
associate M to one of these energy scales. Furthermore, large M values provide an elegant mechanism for generating 
the matter antimatter asymmetry of the Universe 63] . Of course, large M values are experimentally very frustrating. 
It may ultimately prove impossible to experimentally verify whether the New Standard Model is really the correct 
way to describe Nature. 

Here, we explore the opposite end of the M spectrum, M < 1 keV. Such values are phenomenologically very 
intriguing, given that small M values imply the existence of light sterile neutrinos that mix significantly with the active 
neutrinos and can potentially be directly observed. Furthermore, there are several experimental and astrophysical 
phenomena that are best understood if one postulates the existence of light, moderately mixed sterile neutrino. 
We find that by requiring all three right-handed neutrino masses to be less than a few keV we can simultaneously 
explain all neutrino oscillation data, including those from LSND, explain the large peculiar velocities of pulsars, and 
accommodate the production of heavy elements in supernova environments. Our fit also provides constraints for the 
"active" neutrino oscillation parameters, most strongly to the lightest active neutrino mass. All successful parameter 
choices that accommodate the LSND data require m; to be "large" (mi > 0.1 eV), and the "best fit" requires all 
active neutrino masses to be quasi-degenerate. 

On the negative side, low M values are, theoretically, rather puzzling. In order to obtain the observed light 
neutrino masses, neutrino Yukawa couplings are required to be much smaller than the electron Yukawa coupling, and 
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it is tempting to believe that such small numbers are proof that a more satisfying understanding of fermion masses 
must exist. Furthermore, thermal leptogenesis is no longer an option (see, however, [64|]). Finally, the fact that M is 
naively unrelated to other mass scales can also be perceived as disheartening, but in our opinion, should be interpreted 
as evidence that there is more to the lepton sector than meets the eye. 

Regardless of one's preference for a high or low seesaw energy scale, and independent of whether the data from 
LSND and the astronomical observables discussed above have anything to do with sterile neutrinos, our main point 
is that the determination of M is an experimental issue. In the near/intermediate future, low energy seesaw scales 
will be probed by several experiments; most importantly measurements of the end-point of tritium beta-decay, the 
MiniBooNE experiment, searches for neutrinoless double-beta decay and, if we get lucky, the detection of neutrinos 
from a nearby supernova explosion. We find, for example, that Katrin should be sensitive to seesaw energy scales below 
tens of keV if all right-handed neutrino masses are similar, while null results from MiniBooNE would severely constrain 
right-handed neutrino masses below several eV. We conclude by pointing out that larger (but still "small") values 
of M are much harder to constrain. For GeV sterile neutrinos, typical active-sterile mixing angles are < 10~ 10 , 
probably too small to observe in particle physics processes. It is frustrating (and, we hope, ultimately false) that we 
seem to be unable to experimentally distinguish M ~ 1 GeV from M ~ 10 14 GeV. . . 
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